Damping of the de Haas-van Alphen oscillations in the superconducting state of MgB 2 
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The de Haas- van Alphen (dHvA) signal arising from orbits on the n Fermi surface sheet of the two- 
gap superconductor MgB2 has been observed in the vortex state below H C 2. An extra attenuation 
of the dHvA signal, beyond those effects described in the conventional Lifshitz-Kosevich expression, 
is seen due to the opening of the superconducting gap. Our data show that the tt band gap is 
still present up to H C 2- The data are compared to current theories of dHvA oscillations in the 
superconducting state which allow us to extract estimates for the evolution of the it band gap 
with magnetic field. Contrary to results for other materials, we find that the most recent theories 
dramatically underestimate the damping in MgB2. 



I. INTRODUCTION 

The existence of two distinct energy gaps in the su- 
perconducting state of MgB 2 has been demonstrated by 
a number of experiments, including: tunneling, 1 specific 
heat, 2 ' 3 magnetic penetration depth 4 and angle resolved 
photoemission spectroscopy (ARPES). 5,6 Although there 
have been indications of two, or multi-gap, effects in 
other materials, 7,8 it is in MgB 2 where the behavior 
has been most thoroughly investigated. Band structure 
calculations 9 have shown that the Fermi surface of MgB2 
consists of four sheets, mostly arising from the boron or- 
bitals. Two quasi-two-dimensional sheets originate from 
the boron a orbitals and two more isotropic honeycomb 
shaped sheets from the boron 7r orbitals. Theoretical 
work 10,11 has predicted that the superconducting energy 
gap is substantially larger on the a sheets than on the 7r 
sheets (at zero field and temperature the two gaps have 
been measured to be A a ~ 78 K and A w ~ 29 K). 4 ' 12 

Although there have been several studies which have 
accurately measured the temperature dependence of the 
two gaps, 1,13 ' 14 there has been much less work in es- 
tablishing how they evolve with magnetic field. Theo- 
retical work 15 predicts that at high field both gaps de- 
crease towards zero at a common upper critical field 
(H C 2) value, although at lower field is depressed much 
more rapidly than A a , particularly for _ff||a6. 16 Experi- 
mentally, the only direct studies of this have been by 
point contact spectroscopy. It has been suggested 17 that 
A w goes to zero at ~ 1 T (for H\\c), which is much lower 
than H C 2, although others 18,19 have concluded that A„. 
remains finite above 1 T but becomes increasingly diffi- 
cult to resolve because of scattering. In this paper, we 
use dHvA measurements as a probe of the gap in high 
magnetic fields. Our data clearly show the presence of a 
gap on one 7r-band sheet right up to the 'bulk' H C 2- 

The existence of dHvA oscillations in the supercon- 
ducting state has been the subject of study for a number 
of years, and has been observed in a number of different 
materials: NbSe 2 (Ref. 20), V 3 Si (Ref. 21), Nb 3 Sn (Ref. 



22), k-(ET) 2 Cu(NCS) 2 (Ref. 23), YNi 2 B 2 C (Ref. 24) 
and several heavy fermion compounds. 25-29 The effect is 
usually thought to arise through the overlap of quasipar- 
ticle states outside the vortex core. The amplitude of 
the oscillations is governed by the magnitude of the field 
dependent superconducting energy gap. In theory the 
effect can be used to resolve the gaps on different Fermi 
surface sheets and to probe the field dependence of these 
gaps. However, in the current study scattering restricts 
our study to only one Fermi surface sheet. 

The amplitude of the dHvA oscillations is interpreted 
using the Lifshitz-Kosevich equation for the oscillatory 
torque T of a three-dimensional Fermi liquid. The am- 
plitude of the first harmonic is given by 30,31 
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Here the dHvA frequency F is related to the extremal 
area (A) of the orbit in fc-space by F = (fi,/27re)A, 
A" = d 2 A/dk 2 is the curvature factor and ip is the phase. 

The factor R T = A/(sinhA) where X = —^ s - w ^-, 
accounts for the effects of thermal broadening of the 
Landau-levels. It is from this temperature dependent 
term that the quasi-particle effective mass m* is de- 
termined. The Dingle factor accounts for the effect of 
impurities and is given by Rd = exp(— it/uj c t), where 
lo c = eB/rriB, tub is the unenhanced band mass 30 ' 31 and 
r is the scattering time. The spin splitting factor Rs ac- 
counts for the reduction of amplitude caused by beating 
between the spin-up and spin-down Fermi surfaces. It 
is given by Rs = cos(irngrriB(l + S)/2m e ) where S is 
the orbitally averaged exchange-correlation (Stoner) en- 
hancement factor, g is the electron <?-factor, m e is the 
free-electron mass and n is an integer. The final fac- 
tor, RsCi parameterizes the effect of superconductivity 
on the dHvA amplitude. It is this factor which is the 
primary focus of the current work and will be described 
more fully below. 

The calculated Fermi surface of MgB 2 , with 
predicted 32-34 dHvA orbits, is shown in Fig. 1. We de- 
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FIG. 1: Fermi surface of MgB2 with predicted dHvA orbits 
(for frequencies less than lOkT). The orbit which survives into 
the superconducting state, and is the subject of the majority 
of the present study, is orbit 3. [Figure adapted from Kortus 
et al. 9 ] 

note the frequency of each orbit, n, at a general angle by 
F n . Previous studies 35,36 have succeeded in observing or- 
bits 1 to 6, thus verifying the topology of the calculated 
Fermi surface. The measured /c-space areas of the or- 
bits and the quasiparticle effective masses were found to 
be in good overall agreement with the calculations. 32-34 
Orbit 7 (Ref. 32) was not observed experimentally, per- 
haps because of a slight departure of Fermi surface topol- 
ogy from the calculations, or alternatively because of in- 
creased scattering on this orbit. 
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FIG. 2: The main figure shows the amplitude of the oscil- 
lations for orbit 3 at 16. 5T (symbols) as a function of angle 
(9) (the solid line is a guide to the eye). The dashed line 
shows the extrapolated amplitude at H = H C 2(0) using the 
measured Dingle factor Rd- The inset shows the correspond- 
ing angular dependence of the amplitude of orbits 1 and 2 at 
16.5 T fitted according to Eq. 1 (solid lines), and the extrap- 
olation to H = H C 2{9) (dashed lines). The noise level of our 
measurement circuit is indicated on both parts of the figure. 
The mean-free-paths were measured to be 500A, 610A and 
840A for Fi, F2 and F3 respectively. 



The relative amplitude of the various orbits is some- 
what sample dependent. We have found that the mean- 
free-paths for different orbits does not change in a uni- 
form way between samples. 36 Orbits 1-3 have by far the 
largest amplitude for fields below 20 T. The frequencies 
of these orbits all vary approximately like 1 / cos 9 (or 
l/sin#), although for F\ and F 2 there are some depar- 
tures from this simple behavior due to the warping of the 
a band sheet (see Ref. 35 for details). The amplitude of 
orbits 1 - 3 as a function of angle 9 as the magnetic field 
is rotated from B\\c (9=0) to £?||a (0=90°) is shown in 
Fig. 2. The strong angular dependence of H C 2, combined 
with the exponential attenuation of the signal by the Din- 
gle factor, limits the range of angles where we are able to 
observe oscillations below H c2 to 66° < 9 <81° for the 
present crystals. The amplitude of the signal from orbits 
1 and 2 at H C 2 (Fig. 2 inset) is many orders of magni- 
tude below our noise level, and so our study is limited 
to orbit 3 on the electron-like tt sheet. We estimate that 
the mean-free-path I would have to be ~ 3 times larger 
on orbit 1 for us to be able to see these oscillations at 
HqH = 4 T (for these crystals fi H c2 ^ 4.3 T at 9 = 30°, 
although this does decrease with increasing t [Ref. 36]). 



II. EXPERIMENTAL DETAILS 

Rather than measure the oscillations in the mag- 
netic susceptibility, as in a conventional field modulated 
dHvA experiment, we have used miniature piezoresistive 
cantilevers 37 to measure the oscillations in the torque. 
We have found that this latter technique is significantly 
more sensitive than the former for the small size of the 
available high quality crystals of MgB 2 . 

The MgB2 crystals were fixed with epoxy to the end 
of a boron-doped silicon cantilever about 150/im long. 
Deflections in the cantilever were measured through the 
strain-induced changes in its electrical resistance, using 
an AC bridge technique. The torque values are reported 
here in units of bridge resistance R, i.e, the off-balance 
voltage divided by the excitation current (the change 
in cantilever resistance is 4R). We estimate, using the 
weight of the sample, that T ~ 10~ n R (r in Nm and R 
in f2). The noise limit is around 2mfl or ~ 1CP 14 Nm 
(this can vary by up to a factor of 10 between different 
levers) . In addition to the changes in resistance caused by 
the torque there is a small monotonic magnetoresistance 
of the sensor which we have not attempted to correct for. 

The cantilever is mounted on a single axis rotation 
stage in a 3 He cryostat inside the bore of a 19 T supercon- 
ducting magnet (20.5 T at T=2.2 K). Unless otherwise 
stated, all measurements in this paper were performed 
in liquid 3 He at 320±20 mK. The orientation of the can- 
tilever mount is detected using a pickup coil and a small 
modulated field collinear with the DC field. The small 
offset between the crystal plane and the cantilever mount 
was corrected using the symmetry of the dHvA frequen- 
cies around the 9 = 0° and 90° points giving an uncer- 
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tainty in the out of plane angle of ±0.2° (at 9 ~ 90°). The 
in-plane orientation of the crystal, <j>, (measured from the 
a-axis) is fixed for each run and measured by optical pho- 
tographs and Laue X-ray diffraction (to ±5°). A more 
precise determination of <f> was made by comparing the 
minimum in the frequency of orbit 3 as a function of 9 
to the known frequency for this orbit at 9 — <p = 0°, 
as measured on a large number of other crystals 38 (the 
minimum frequency varies approximately as 1 / cos <fi) . 

The torque, both due to the oscillations and the back- 
ground magnetism, causes a deflection of the lever, and 
so 9 is not quite constant during the field sweep. Us- 
ing the strong angular dependence of the dHvA frequen- 
cies, we are able to measure the correspondence between 
the measured torque signal and the angular deflection of 
the cantilever. We find that the deflection is given by 
A9 = (0.05 ± 0.01)°/i7, which is approximately 4 times 
smaller than the value estimated from the manufacturer's 
data sheet. 37 This 'torque interaction' effect gives rises to 
spurious harmonics and mixing of the main frequencies. 
As the amplitude of the signal is also angle dependent 
(see Fig. 2) the field dependence of the signal will not ex- 
actly follow Eq. 1. Fortunately, this effect is small enough 
to be neglected for the purposes of the current work. This 
was verified by measuring samples with different masses 
and comparing up and down field sweeps (see later). 

The samples used in this work were grown by a high 
pressure synthesis route as described in Ref. 39. Most 
of the work reported here was conducted on a crystal 
measuring 300 x 160 x 30^m 3 (mass=3.9/ig). The T c 
of samples from this batch (AN77) was measured to be 
38.5 K. 



III. RESULTS 

In Fig. 3 we show the measured torque versus field 
over a large range of field for a small sample from the 
same batch as our main crystal (mass=0.4/xg). This small 
sample was selected for this part of the study to avoid 
overstressing the cantilever. The general shape of the 
curve, a bell shaped bump with peak at around H c2 /2, 
has been explained theoretically in Ref. 40. Near H C 2 
there is evidently a pronounced peak effect which grows 
in size as the angle is increased. In addition, there is a 
significant region above the peak effect region where the 
torque is sizeable and remains hysteretic. Similar fea- 
tures have been reported previously by Angst et al. 41 ' 42 
in their torque study of MgB 2 single crystals. Their study 
was conducted at higher temperature and lower field than 
the current work, and the peak effect they observed was 
less marked. In the inset to Fig. 3 we show data for the 
ac susceptibility of another crystal from the same batch 
(mass=17/ig). Many of the same features evident in the 
torque data are also visible here, namely the broad super- 
conducting transition and the pronounced peak effect. 

The existence of a pronounced curvature in the region 
of the normal/superconducting transition along with the 
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FIG. 3: Torque versus field at several values of at T 
=320 mK. Both up and down field sweeps are shown. Inset: 
AC susceptibility (real part) versus field for different crystals 
from the same batch (at T—1.2 K). The excitation field was 
38 Oe at a frequency of 72 Hz. 



peak effect feature means that there is no unambiguous 
way of extracting H c2 from the measured torque curves. 
The rounding of the transition is common in high T c ma- 
terials where it is usually attributed to the presence of 
thermal or quantum fluctuations, which are enhanced rel- 
ative to conventional low T c materials due to the low di- 
mensionality and high T c of the former. MgB2 is not 
strongly anisotropic, but has a relatively high H C 2 (for 
H parallel to the basal plane) and hence small coherence 
lengths. The existence of hysteresis in this rounding re- 
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FIG. 4: H* 2 as determined by torque measurements and AC 
susceptibility measured at 320mK (solid symbols) , along with 
a fit to Eq. (2). The open symbols show the values of H* 2 (9) 
extracted from the fits of Rsc using Eq. (3) for two different 
values of <b and are discussed later. 
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FIG. 5: Raw torque versus field close to H C 2 for several angles 



gion however, points to an explanation either in terms 
of surface superconductivity 43 or an as yet unexplained 
two gap effect, rather than fluctuations. It may also be 
possible that it could result from crystal inhomogeneity, 
although this is unlikely as reproducible behavior is found 
between crystals from the same batch. 

We find that the angular dependence of H c2 derived 
by various extrapolation schemes (e.g., position of peak 
effect maximum or a linear extrapolation from low field) 
gives different absolute values, but essentially the same 
form for H c2 (8)/H c2 (9 = 0) (the same conclusion was 
found by Angst et al. 42 ). As we shall see later, the onset 
of damping of the dHvA oscillations essentially coincides 
with a critical field H * 2 derived by extrapolating the low 
field behavior, and so it seems likely that this represents 
the bulk upper critical field, at least for the n band. 

In Fig. 4 we show values of H* 2 extracted from the 
torque data in Fig. 3. These values are in good agree- 
ment with those extrapolated from higher temperature 
studies. 41 In this figure we also show H* 2 extracted from 
ac susceptibility data (inset to Fig. 3) in a similar way. 

In an anisotropic superconductor the angular depen- 
dence of H c2 is usually described by 



H c2 (9) 
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where jh is the anisotropy of H c2 . In MgB 2 , the contri- 
bution of multiple Fermi surface sheets with different su- 
perconducting gaps is known to cause jh to increase with 
decreasing temperature and to cause deviations from the 
angular dependence predicted Eq. (2). 42 However, these 
are small on the scale of Fig. 4, and are most pronounced 
near 9 = 90°. A fit of this to the data is shown as a solid 
line in Fig. 4 (j H = 7 ± 0.5 and = 3.7± 0.3 T ). 

Torque versus field for our main crystal is shown in 
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FIG. 6: Amplitude of dHvA oscillations divided by Rt 
versus inverse field, at several different values of 9 [<f> — 5°]. 
The data have been offset for clarity. The actual variation of 
amplitude with angle is shown in Fig. 2. 



Fig. 5 for three angles (9 = 63.8°, 69.7° and 72.5°). 
At these angles only one dHvA frequency (F 3 ) is visi- 
ble in our field range. For 9 = 69.7° pronounced hys- 
teresis is observable which disappears at an irreversibil- 
ity field /io-ffi rr ~ f5.5 T. There is a small residual 
field dependence in the background torque which ex- 
tends to ~16 T. We note that both these fields are far 
in excess of our estimate of H c2 obtained from the raw 
torque {^oH* 2 ~ If T). Similar behavior can be seen for 
9 = 72.5°. 

The field dependent amplitude A of the dHvA oscilla- 
tions was extracted by fitting A sin(27rF 3 /B + ip) + aB + b 
(the linear term accounts for the background torque) to 
different sections of data comprising of l\ oscillations. 
The data were then divided by the weakly field dependent 
term B^Rt to give A oc RdRsc (see Eq. 1). The quasi- 
particle effective mass m* in the expression for Rt was 
determined by measuring the temperature dependence of 
the dHvA amplitude [for F 3 , m* = (0.456 ± 0.005)m e at 
9 = 70.8°]. 

In Fig. 6 we show A versus inverse field on semi-log axes 
('Dingle plot') for several angles. For 9 = 63.8°, A varies 
strictly exponentially with inverse field, A oc exp(— 74/_B) 
from which we estimate that the quasiparticle mean free 
path on this orbit is 840 ± 20 A. For this angle there 
is no evidence of superconductivity in cither the back- 
ground or oscillatory torque (see Fig. 5). As 9 is increased 
H c2 increases sharply and the Dingle plots have marked 
downturns below some critical field. We attribute this 
to the opening of the gap as the sample enters the su- 
perconducting state. To extract Rsc wc fit A(l/B) in 
the normal state to the exponential expression for Rd 
and extrapolate this dependence into the superconduct- 
ing state (solid line in Fig. 6). Dividing A by Rd thus 
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FIG. 7: Raw torque and superconducting damping, Rsc at 
9 = 69.7°, for both up and down field sweep directions (indi- 
cated by arrows). 



yields Rsc- An example of this is shown in Fig. 7 in 
which we also show the raw hysterctic torque. From this 
figure it can be seen that the onset of attenuation of the 
dHvA signal from R§c does not occur until the field is 
below .ffi rr , and closely corresponds to the H* 2 value (in- 
dicated by an arrow on Fig. 7) deduced by extrapolating 
the lower field data. 

Most previous studies of the dHvA effect in the super- 
conducting state have been conducted by the field mod- 
ulation technique. In studies where this technique has 
been used the strong flux pinning close to H C 2 (peak ef- 
fect) has prevented data collection in this region. This 
did not present a serious problem as the oscillations 
could be observed far below iJ c2 [e.g., down to -ff c2 /5 in 
YNi 2 B 2 C (Ref. 44) or H c2 /2 in V 3 Si (Ref. 45)]. Studies 
of YNi 2 B 2 C have shown that, in contrast to field mod- 
ulation measurements, torque measurements 46 are not 
affected by this increased pinning. In MgB 2 , oscillations 
are only observable very close to H c2 and so all our mea- 
surements are essentially carried out in the peak effect 
region. 

For <f> near to 0° there is a pronounced dip is found in 
the amplitude of the dHvA signal at 9 ~ 75° (see Fig. 
8). This feature has been present in every crystal we 
have measured to date and always occurs when 9 and (f> 
are adjusted so that F 3 ~ 2800 T. Previously 35 ' 36 this 
feature has been attributed to a 'spin-zero', i.e., an an- 
gle where the angle dependent band mass multiplied by 
the Stoner factor exactly equals m e /2, so that Rs = 0. 
However, a detailed analysis of the field dependence of 
the dHvA amplitude for angles above 9 ~ 70° reveals 
that in fact this dip is produced by a beat between two 
dHvA orbits with very similar frequencies. The inset to 
Fig. 8 shows Dingle plots close to 8 ~ 75°. The data 
clearly show oscillatory beating damped by the Dingle 
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FIG. 8: Amplitude versus 9 at (M)H — 16.5 T, for two different 
in- plane rotation angles = 5° and 24°. Open symbols are 
the amplitude at 18. 5T scaled by the appropriate Dingle fac- 
tor. Inset: Dingle plots for 9 near the minimum in amplitude 
at 75° (Data have been offset for clarity) . 



factor. A fit that takes account of this effect is shown for 
9 = 75.1° in Fig. 8. The frequency difference is - 27 T at 
9 = 90° and extrapolates to zero at 6 = 70° (for 0-0°). 
The relative amplitudes of the two frequencies are almost 
equal over the whole angular range. The characteristics 
of this orbit are very similar to an additional orbit pre- 
dicted by Harima. 32 A full account of this result will be 
given elsewhere. 47 

A consequence of the above dip feature is that at angles 
close to or beyond 9 ~ 75° it is difficult to extract Rsc 
as the functional form of A(l/B) is more complicated. 
For this reason we performed a second set of sweeps at 
a different in- plane angle 4>. With cj> = 24°, F 3 (9 = 0) is 
increased to 2930 T and thus the dip does not occur for 
any value of 9 (see Fig. 8). In Fig. 9 we show the A(l/B) 
curves for this in-plane angle. 

As mentioned earlier, the deflection of the cantilever, 
either by the background or oscillatory torque, means 
that the fields sweeps are not done at strictly constant 
angle. As the dHvA amplitude is angle dependent this 
may cause some additional field dependence to A, which 
becomes large as the sample enters the superconducting 
state, and might cause error in our determination of Rsc- 
The insensitivity of our results for Rsc to this deflec- 
tion is perhaps best demonstrated by comparing A(l/B) 
curves for up and down field sweeps. In Fig. 10 we show 
both the raw torque signal and A(l/B) for 9 = 71.7°. 
There is a large difference in the background torque be- 
low H 1TT in the two cases, which means that the sample 
is deflected in opposite directions (maximum deflection 
~ 0.7°), however it can be see that the drop in A(l/B) 
below H c2 is virtually identical (see also Fig. 7). As a fur- 
ther check we have repeated some our measurements on a 
fragment of crystal cut from our main sample which had 
mass of only 1.3^g (~ 0.3x mass of main sample). We 
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FIG. 9: Dingle plots for in-plane rotation cj> = 24°. The solid 
lines are fits to the Dingle factor Rd [Eq. (1)] for the data 
above H C 2- 



found identical behavior, showing again that the lever 
deflection does not affect our results. The close corre- 
spondence of the results for Rsc for up and down sweeps 
is also a strong indication that our data is not strongly 
affected by pinning in the peak effect region. 



IV. EXTRACTING THE GAP FROM THE 
SUPERCONDUCTING STATE DAMPING 
FACTOR 

dHvA oscillations in the superconducting state have 
been observed in many different materials, and in all 
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FIG. 10: Raw torque and A versus inverse field at 6 — 71.7° 
for both up and down field sweeps. The direction of the field 
sweep is indicated with the arrows. 



cases the oscillations persist into the superconducting 
state with the same frequency as in the normal state 
but with reduced amplitude. As the dHvA amplitude 
depends exponentially on field, all the materials inves- 
tigated to date are those with high H c2 values, and are 
not necessarily conventional. For example, NbSe2 and 
Nb 3 Sn, have recently been suspected of having an exotic 
gap structure. 48,49 In some ways, MgB 2 is the best un- 
derstood of all the materials where this phenomena has 
been observed to date, although its multiple energy gap 
structure may complicate the analysis. 

Several theoretical models have been proposed to de- 
scribe the effect (see Ref. 50 and references therein). 
Two theories have proved most successful in describing 
the data to date. Maki 51 used a semiclassical approach 
following that of Brandt et al. 52 , in which the gap is 
approximated by the spatially averaged value of A 2 and 
the magnetic field is considered to be uniform. For quasi- 
particles moving perpendicular to the magnetic field this 
model predicts that the excitation spectrum is gapless. 
ft is this gaplessness which is the physical origin of the 
quantum oscillations. The extra damping in the super- 
conducting state is given by 



R sc = cxp 



Ag(g) 
frw c 



(3) 



In this expression, the effective field dependent energy 
gap Ae{B) is resolved to that on a particular quasipar- 
ticle orbit. An equivalent expression was also derived 
by different authors starting from alternative physical 
pictures. 31,53 

An alternative expression was derived by Miyake, 54 
who considered the effect of incorporating the zero field 
BCS quasiparticle energy and occupation numbers into 
the usual Lifshitz-Kosevich theory. In the superconduct- 
ing state the sharp step in the Fermi function is replaced 
with the BCS quasiparticle occupation function (|ufe| 2 ) 
whose width is set by the superconducting energy gap. 
This gives a more gradual emptying of Landau levels than 
in the normal state, and hence a reduced dHvA ampli- 
tude. The size of the damping factor in this model is 
given by 



Rsc — xK\{x), x = 2tt' 



hu> r 



(4) 



where K\{x) is the Bcsscl function of the second kind. 
Miller and Gyorffy 55 have derived the same result start- 
ing from the Bogoliubov-de Gennes equations. 

Studies in other materials have shown that Ae(B = 0) 
extracted by fitting data to these expressions 45 does not 
always coincide with the known energy gap measured by 
more conventional means (for example tunneling). The 
correspondence between and the superconducting 
gap is different in each theory, and this may also be differ- 
ent for different materials. For example, Janssen et al. 45 
found that for NbSe2 the Maki model gave A#/A = 0.63 
whereas the Miyake model gave, A#/A = 0.11. For 
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TABLE I: Parameters extracted from the fits to Rsc 
Theory <f> A E (B = 0) fO 3 x a 



Maki 


5° 


200±40 K 


55±f0 


Maki 


24° 


320±40 K 


34±f0 


Miyake 


5° 


32±f0 K 


27±5 


Miyake 


24° 


60±f0 K 


22±5 



VsSi, Ae/A was found to be 2.6 and 0.61 for the two 
models respectively. It is generally found that the Maki 
model gives values for A E which are much larger than 
the Miyake model. 

To compare these theories with experiment we can take 
two approaches. The first is to assume a form for the field 
dependence of the effective gap and then fit Eq. (3) or Eq. 
(4) to Rsc(B) directly. The second is to solve Eq. (3) or 
Eq. (4) for Ae{B) at each field point and then compare 
Ae(B) to the expected behavior. We shall show both 
methods below. 

In the usual mean-field BCS theory we expect 



A|(g) 
A|(0) 



= 1 



B 

B C 2 



(5) 



However, the observed rounding of the superconducting 
transition means that dA/dH does not change discon- 
tinuously at H c2 and so a more appropriate expression 
is 



A|(B) 
A|(0) 



B 

B C 2 



+ 1 



This form of A(B) was used by Clayton et al 23 to ac- 
count for the effect of superconducting fluctuations, as it 
interpolating the mean field result at low field to the form 
expected for a strongly fluctuating system near H c2 . The 
phcnomcnological parameter a sets the strength of the 
fluctuations. As it is not clear where the rounding of the 
superconducting transition in MgB 2 is caused by fluctua- 
tions, we will regard a as a phcnomcnological parameter 
which describes the rounding, rather than attributing to 
it for now any particular physical significance. In any 
case, the values of Ae(B = 0) are mostly determined by 
the portion of the Rsc curve well below H c2 (i.e., the 
data for the higher values of 9). 

Combining Eq. 6 with either Eq. (3) or Eq. (4) results 
in 3 fitting parameters per curve (A^(0), H c2 and a). 
We found that in general there was too much covariancc 
between the parameters to arrive at accurate values for 
Ae(0) if we allowed all three to vary between the fits 
at each angle. As we do not expect either Ae(0) or a 
to vary strongly with 9 we found that more consistent 
results were obtained by fitting to the data for Rsc(B) 
at all angles simultaneously, only allowing H c2 to vary as 
a function of angle. 




11 13 
H H [T] 



FIG. 11: Rsc for <t> = 5°(top panel) and <j) = 24° (bottom 
panel). The solid and dashed lines are fits to the Maki (Eq. 
3) and Miyaki (Eq. 4) theories for the superconducting state 
damping 



The fits obtained for both <fr = 5° and <f> = 24° are 
shown in Fig. 11. In these fits we have taken mj(9) = 
0.3l5m e F 3 /F°, i.e., assuming the band mass scales like 
the dHvA frequency. 34 - 56 The values of A E (B = 0) and 
a and are shown in Table I. It can be seen that both the- 
ories provide an adequate fit to the experimental data. 
The main difference is in the size of the extracted su- 
perconducting gap. As orbit 3 is on the electron like tt 
sheet of Fermi surface we should compare these values 
with the zero field value for the smaller MgB 2 gap (A s ~ 
29 K). 4 The Miyake fit gives values comparable to this 
whereas the Maki fit gives values up to 10 times larger. 
The trend is similar to that found for NbSe2 . The values 
of a are also somewhat larger in the Maki fits to those 
in the Miyake fits. The values of H c2 extracted from our 
fits are shown in Fig. 4. The values are very similar to 
those extracted from the background torque, consistent 
with our interpretation of our extrapolated H* 2 as the 
bulk critical field for the tt band. There is a consistent 
difference in H c2 (9) values for the two values of <f>, with 
the values for cf> = 5° being <~ 10% higher than those for 
<t> = 24°. This implies a slight in- plane anisotropy of H c2 . 

Although there is some ambiguity between the theories 
as to the magnitude of the gap, it is clear that a sizeable 
gap does exist on this sheet of Fermi surface right up to 
H c2 . There is some difference between the gap values de- 
rived from the data at the two in-plane rotation angles. 
We note that although there is a large difference between 
the gaps on the a and 7r Fermi surface sheets, the differ- 
ence within a sheet (or indeed pair of sheets) is expected 
to be small. 57 One reason for this orientation dependence 
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H H[T] 



FIG. 12: Calculation of the field dependence of the supercon- 
ducting gap, by inverting the Rsc curves in Fig. 11 using the 
Maki [Eq. (3)] model. The solid lines are fits to the field de- 
pendent gap [Eq. (6)] (see text). The inset shows the variation 
of the fit parameter a with 8. 



is that the data for </> = 24° cover a wider range of Hj H c2 
and hence, the extrapolated gap value is more accurate. 
For the runs with different <\> where the H c2 values are 
close, the Rsc curves are very similar, and any difference 
results from the field range of the fit. Next, we show in 
Fig. 12 the field dependent effective gap calculated by 
inverting Eq. (3), along with fits to the field dependent 
gap given by Eq. (6). The main difference between using 
Eq. (4) rather than Eq. (3) is just in the magnitude of 
the extracted Ag rather than its field dependence. In 
the fits we have fixed Ae(B = 0) to 300 K for all values 
of 9 and and allowed a to vary. The curvature near 
H r2 and the approximately linear behavior of A 2 E (B) at 
lower field are evident. The covariance between the gap 
and H c2 mentioned above arises from the rounding at 
the transition and the small field range available at some 
angles. There is however a discernible progressive change 
in shape of the Ae(B) curves as 9 and <j) are varied. This 
is reflected in the angular dependence of a in the inset to 
the figure. This may be an artifact of the data analysis 
but could also result from a two gap effect. 

It is clear that the extrapolation to the zero field gap 
is not a trivial one as our data does not extend below 
0.85-ff C 2- This is particularly significant in a material like 
MgB2, where non mean-field gap behavior is expected. 
Numerical calculations 15 ' 16 of the field dependence of the 
energy gaps in MgB2 show that they differ considerably 
from the mean field behavior given by Eq. (5). Close to 
H c 2, for iJ||a&, we find that the numerical data can be ap- 
proximated by Al(B) = 0.24A2(0)(1 - H/H c2 ) 1A . This 
accounts naturally for some (but not all) of the round- 
ing in our Rsc plots. The steep drop of A 2 at low field 
found in the calculations means that a linear extrapo- 
lation from high field underestimates the true zero field 
gap by approximately 50%. 



Recently, there has been further theoretical work 
on the dHvA effect in superconductors. Duncan and 
Gyorffy 58 have extended the work in Rcf. 55 and given 
a new formula for Rsc- We have fitted our data to this 
formula and find fits that are comparable in quality to 
those in Fig. 11, with A E (B = 0) = 200 ± 20 K, for 
(f) = 5°. This is very close to the values calculated from 
the Maki formula [Eq. (3)]. Yasui and Kita 50 have made 
a detailed numerical investigation of the approximations 
used in many of the other previous superconducting state 
dHvA theories. They proposed an equation for Rsc 
which should give quantitative values for the energy gap 
which can be compared to the actual thermodynamic 
values. Their expression for Rsc is 



Rsc = exp 



-27T/3 



A E (B) 
huj c 



(7) 



where (3 — 0.0625 is a numerically evaluated constant. 
This was shown to correctly describe the damping in 
Nb 3 Sn and NbSe2 with A# approximately equal to the 
actual superconducting energy gap. It can be seen that 
this is almost exactly the same as Maki's formula, Eq. (3), 
except for a factor 2(3 /{ir 2 (F/B)?] which equals 0.94 for 
F = 2700 T and B = 15 T. Hence, the gap values ex- 
tracted by fitting this formula are virtually identical to 
those from Eq. (3) i.e., up to 10 times larger then the gap 
values extracted by other means. 



V. CONCLUSIONS 

In conclusion, we have made a detailed study of the 
attenuation of the dHvA signal in MgB 2 as it enters 
the superconducting state. Only a single orbit on the 
electron-like ir band is observable. The data clearly show 
that a sizeable gap exists on this orbit even at high field. 
The transition from normal to superconducting states as 
seen in both the background magnetisation as well as the 
dHvA amplitude are rather broad, which may result from 
either fluctuation effects or a two-gap effect. The data 
provide a new test for theories of the dHvA effect in the 
superconducting state being applied for the first time to 
a material which clearly has two superconducting gaps. 

Both the Miyake and the Maki theories fit the data 
well, although the Miyake theory produces gap values 
which are much closer to those expected from the low 
field data for the tt band. The most recent theoreti- 
cal work 50 ' 58 however, suggests that the Miyake theory 
should seriously overestimate the damping and that the 
Maki model should provide an accurate quantitative esti- 
mate of the average gap on the dHvA orbit. Our data is 
in serious disagreement with this as the gap values found 
from the Maki model fits arc more than a factor 10 times 
larger than the expected zero field tt band gap. We con- 
clude that the present theories of the dHvA effect in the 
superconducting state, when applied to MgB2, dramati- 
cally underestimate the damping. 
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